Gravitational Radiation from Primordial Helical MHD Turbulence 
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We consider gravitational waves (GWs) generated by primordial inverse-cascade helical magneto- 
hydrodynamical (MHD) turbulence produced by bubble collisions at the electroweak phase tran- 
sitions (EWPT). Compared to the unmagnetized EWPT case, the spectrum of MHD-turbulence- 
generated GWs peaks at lower frequency with larger amplitude and can be detected by the proposed 
Laser Interferometer Space Antenna (LISA). 

PACS numbers: 98.70.Vc, 98.80.-k 



When detected, primordial cosmological GWs will pro- 
vide a very valuable probe of the very early Universe 
[H|. Various mechanisms that generate such GWs have 
been discussed: quantum fluctuations 0; bubble wall 
motion and collisions during phase transitions [A ]M ; cos- 
mic strings cosmological magnetic fields II, Mil; and 
plasma turbulence 



11, |12H. From the direct de- 



nologic 



tection point of view, GWs generated during the EWPT 
are promi sing since their peak frequency lies in or near 
the LISA [13( frequency band 14j, however, to produce 
a detectable signal the EWPT must be strong enough 



Il7j |. Currently discussed EWPT models do not 
predict an observable GW signal from bubble collisions 
17 1 , nor for GWs produced by unmagnetized turbulence 



Here we study the generation of GWs during a first- 
order EWPT assuming that bubble collisions produce he- 
lical MHD turbulence. 1 In the case of unmagnetized hy- 
drodynamical turbulence the peak frequency of the GW 
power spectrum is determined by the inverse turn-over 
time of the largest eddy and the energy-scale when the 
GW is generated. Recently discussed modifications of the 
standard EWPT model place the transition at a higher 



energy-scale (16J. As a result, the GW power spectrum 



peak frequency is shifted to higher frequency which, since 
the GW spectrum is sharply peaked, reduces the possibil- 
ity of detection by LISA. On the other hand, in the case 
of MHD turbulence the presence of an energy inverse- 
cascade leads to an increase in the effective size of the 
largest eddy (now associated with an helical magnetic 



1 Kinetic or magnetic helicity generation at the EWPT is studied 
in Refs. [l8| . Previously we studied generation of GWs by direct- 
cascade turbulence and found that, due to parity violation in 
the early Universe, the induced GWs are circularly polarized 
[l ll | . Polarized GWs are present in other models [l9| . and the 
polarization of the GW background is in principle observable, 
either directly or through the CMB fHw\. 



field), and can result in the GW power spectrum peak- 
ing in the LISA band, with amplitude large enough to be 
detected by LISA. We adapt the technique developed in 
Ref. [13] to study this case here. We model MHD turbu- 
lence and obtain the GW spectrum by using an analogy 
with the theory of sound wave production by hydrody- 



namical turbulence 

Since the turbulent fluctuations are stochastic, so are 
the generated GWs. The GW energy density is [141 ] 
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Here the times t'^ = t — |x — x'(")|, i and j are spa- 
tial indices (repeated indices are summed), the source 
Sij(x,t) — Tij(x,t) — SijT^(x,t)/3 is the traceless part 
of the stress-energy tensor T%j , G is the gravitational con- 
stant, and we use natural units with h = 1 = c. We as- 
sume that the turbulence exists for a time short enough 
to neglect the cosmological expansion during GW pro- 
duction. We consider metric perturbations in the far- 
field limit (i.e. for x 3> d, where d is a characteristic 
length-scale of the source region), where GWs are the 
only metric perturbations [26j], and replace |x — x'| by |x| 
in Eq. {T]). If the turbulence is stationary then the GW 
spectral energy density /(x, u>) (pgwJx) = J dw/(x, ui) 
where uj is the angular frequency) is [121 ] 
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x , —u>, ,j j . (2) 



Here Hijij(x' , k, w) (where k is a proper wavevector) is 
the (double traced) four-dimensional Fourier transform of 
the two-point time-delayed fourth-order correlation ten- 
sor, (Sij(x.\t)Si m (x.",t + T))/w 2 , with respect to x" — x' 
and r, where w = p+p is the enthalpy density and p and 
p the pressure and energy density of the plasma. 

We assume that primordial MHD turbulence is gener- 
ated at time at proper length-scale — 27r//co with 



2 



characteristic velocity perturbation vq. 2 The dynamics of 
MHD turbulence is dominated by Alfven waves for which 
the magnetic and kinetic energy densities are in approx- 

In this case v 



27|. 



bo, where 



imate equipartition 
60 = £>o / V 47rw is the characteristic magnetic field per- 
turbation expressed in velocity units. While MHD turbu- 
lence is isotropic on large scales, it is locally anisotropic 
on small scales [H} , resulting in small-scale anisotropy in 
the generated GW background. However, GWs are gen- 
erated mainly by the largest eddies 12j so we adopt an 
isotropic turbulence model, and thus the magnetic field 
two-point correlation function is ( 6* f k, t)bj(k' , t + r)) = 
F^(k, t)f( V (k), r)S(k - k'), with O 



E (k, t) 
ink 2 



iSijiki 



H M (k,t) 
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Here P y (k) = 5 lj ~ hk^/k 2 , and E M (k,t) and H M (k,t) 
are the magnetic field energy and helicity densities. The 
Schawarz inequality implies \H M (k,t)\ < 2E M (k,t)/k 
[27j| . For the total magnetic energy £m(*) = / E M (k,t)dk 
and helicity H M (t) = j H M (k,t)dk we get H M {t) < 
2£ M (t)£ M (t), where f M (t) = J E M (k 1 t)k- 1 dk/£ M {t) is 
the magnetic-eddy correlation length. r){k) is an auto- 
correlation function that determines the characteristic 
function f(r)(k),T) that describes the temporal decorre- 
lation of turbulent fluctuations. In the following we use 
/(^(fc),r)=exp(-^ 2 (fc)r 2 /4) HI. 

After generation primordial turbulence freely decays. 
We adopt the decaying MHD turbulence model of Refs. 
M, [HI. For non-zero initial magnetic helicity turbulence 
decay is a two stage process. First decay stage dynamics 
is governed by a direct cascade of energy density lasting 
for a time r s o = sqTq, several times (sq ~ 3 — 5) longer 
then the characteristic largest-eddy turn-over time tq = 
lo/vo — 27r/fcoi>o- During the first stage energy density 
flows from large to small scales and finally dissipates on 
scales ~ Id — 2ir /kd (kd S> ho) where one of the Reynolds 
numbers becomes of order unity. Due to the selective de- 
cay effect [27| magnetic helicity is nearly conserved dur- 
ing this stage [3l| . To compute the GWs generated by de- 
caying MHD turbulence, we assume that decaying turbu- 
lence lasting for time t s q is equivalent to stationary tur- 
bulence lasting for time r s o/2. This can be justified us- 
ing the Proudman [22I [25| argument for (unmagnetized) 
hydrodynamical turbulence. Consequently, when com- 
puting the emitted GWs we ignore the time dependence 
of E M (k,t) and H M (k,t). We also assume small initial 
magnetic helicity, a* = Hm(U)/[2£ m (.U)£ m (U)] < 1. 



2 We assume that the usual and magnetic Reynolds numbers are 
much greater than unity on scales ~ Iq , otherwise there is no tur- 
bulence. Throughout this paper the symbol ~ represents equal- 
ity to the accuracy of a dimensionless multiplicative constant of 
order unity. 



For E M [k,t) and rj{k) we use the Kolmogorov model, 
£ M (fc,i) ~ £ 2/3 fc" 5/3 , V(k) = e l ' 5 k 2 ' 3 /^, (4) 

for ko < k < kd- Here e ~ fco^o ^ s ^ ne energy dissipation 

rate per unit enthalpy. 

At the end of the first stage turbulence relaxes to a 

maximally helical state, a s0 ~ 1 [U, [32|. Accounting 

for conservation of magnetic helicity, the characteristic 

velocity and magnetic field perturbations at this stage are 

1/2 1/2 
Vi ~ aj vq and b\ ~ a J bo. Second stage dynamics is 

governed by a magnetic helicity inverse cascade. If both 

Reynolds numbers are large at the end of the first stage, 

magnetic helicity is conserved during the second stage. 

The magnetic eddy correlation length evolves as £m (*) ~ 

, where t\ 



lo \J\ + tprx [30J, l31| where t\ ~ Iq/vx — to/ \fa* is the 
characteristic energy containing eddy turn-over time at 
the beginning of the second stage. The magnetic £m(£) 
and kinetic £it(t) energy densities evolve as (30l. l3l| 



£ M (t) cx (l+t/nj-^wft?, 
£ K (t) cx (1 + t/r 1 )- 1 wvl 



(•5) 



These imply that the characteristic turn-over (r to ) and 
cascade (r cas ) timescales evolve as 



nil + t/n). 



(6) 



To compute the GWs emitted during the second stage 
we use the stationary turbulence model that has the same 
GW output. Introducing the characteristic wavenum- 
ber fcf(t) = 27t/£m(£) and using Eqs. ((5|) we find £m ~ 
wv 2 k^(t)/k and £k ~ ww 2 [fc^(i)/fc ] 2 since b\ ~ V\. 
The time when turbulence is present on scale £m(£) 
is determined by Eq. ([6]) which can be rewritten as 
Teas ~ Ti [ko/k^(t)] 2 . So instead of considering decay- 
ing turbulence, we consider stationary turbulence with 
a scale-dependent duration time (time during which the 
magnetic energy is present at the scale), t s1 ~ Ti[k /k] 2 
(for k = k^ this coincides with r cas ). 

The expression for £m yields the time-independent 

E M (k,t) = dvf/ko = kH M (k,t)/2, k s <k< k . (7) 

Here C\ is a constant of order unity, ks is the small- 
est wavenumber where the inverse cascade stops, and 
the second equation follows from saturating the causal- 
ity condition. For the second stage autocorrelation func- 
tion, which is inversely proportional to the turn-over time 
(|6|), we assume rj(k) = {k/ko) 2 /^/2ttt\. At the largest 
scales there is no efficient dissipation mechanism, so the 
inverse cascade will be stopped at scale ls{t) = 2tt /kg 
where either the cascade timescale r cas reaches the ex- 
pansion timescale H^ 1 = i7 _1 (t^), or when the char- 
acteristic length scale £m(0 ~ Is reaches the Hubble 
radius. These conditions are a+ 1 ^ 2 l 2 s /vok) < H+ 1 or 
Is < H+ 1 (the cascade time is scale dependent and max- 
imal at k = ks)- Defining 7 = Iq/H^ 1 ( 7 < 1), it is 



easy to see that the first condition is fulfilled first and 
consequently ko/ks < («o/7) 1/,2 ci!* • To have an inverse 
cascade requires ko/ks > 1, leading to a constraint on 

1/2 

initial helicity, 7 < MaJ (where M = vq is the turbu- 
lence Mach number). 

The magnetic field perturbation stress-energy tensor 
is Ty r (x,t) = w&j(x, t)6j(x, i). For the first decay stage 
we compute for this magnetic part and then double the 
result to account for approximate magnetic and kinetic 
energy equipartition for Alfven waves. During the sec- 
ond stage, according to Eqs. kinetic energy can be 
neglected compared to magnetic energy. To compute 
Hjjn (k, uj) we assume Millionshchikov quasi-normality 



[22| and use the convolution theorem (for details see Sec. 
Ill of Ref. 12]). Using the (k — > 0) aero-acoustic approx- 
imation, which is accurate for low Mach number [M < 1, 
(and slightly overestimates GWs amplitude for the Mach 
number approaching unity (M — > 1) [12j), we find 



Hijij{k,u}) w Hijij(0,(d) 



(8) 
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The integral is dominated by the contribution of large 
scale (fc ~ ks) perturbations and is maximal at ~ 
aj 2 Mkg/k = 2-kH*. For the first-stage direct-cascade 
turbulence the peak frequency is u£L ~ fc M G3. To 
determine the peak frequency at the current epoch we 
need to account for the cosmological expansion which 
decreases the GW amplitude and frequency by the factor 
a*/ao, where a* and ao are the values of the cosmological 
scale factor at the GW generation and current epochs 

The total GW energy spectrum at a given space- 
time event is obtained by integrating over all source re- 
gions with a light-like separation from that event, and 
includes contributions from GW generated during the 
first and second stages. For the first stage (with du- 
ration time = SnT ) pgyyiid) is given by Eqs. (21) 
and (A3) of Ref. [12j. For the second stage contribu- 
tion we must account for the scale dependence of the 
cascade time. The total GW fractional energy den- 
sity parameter at the moment of emission £7q\y * is 



(m) 



H^^id^/Hl 0. Here the index 



105ff*w 3 E m r : 

to runs over / and II for the first and second decay 
stages, Hq is an angular frequency at the moment 
of emission. The current GW amplitude is related to 
the current fractional energy density parameter through 

h c (f) = 1.26 x 10- 18 (lHz//) [h 2 n GW (f)] 1/2 (where 
ho is the current Hubble parameter Hq in units of 100 
km sec~ ^pc -1 ) (HI, and 



hcU) 



2 x 10" 



100 GeV 



E 



4 cj-(m) 
ijij 



100 

9* 

(o,< 



1/3 



1/2 



(9) 




FIG. 1: The spectrum of gravitational radiation from MHD 
turbulence for g* = 100, T, = 100 GeV, 7 = 0.01, and 
M — l/\/3, for four different initial magnetic helicity values, 
a* = (solid line), a* = 0.02 (dashed line), a* = 0.05 (dash- 
dotted line), and a* =0.1 (dotted line). The bold line is the 
LISA design sensitivity curve. 



Here the linear frequency / = (a*/ao)/* with /* = 
w^/27r, T+ and <?* are the temperature and effective num- 
ber of relativistic (all fields) degrees of freedom at scale 
factor a*. 

Figure 1 shows hc(f) for a few initial magnetic he- 
licity values. GWs emitted during direct-cascade un- 
magnetized turbulence peak at current /mi x ~ Mv* 
[lij ]. We find that the MHD-inverse-cascade gener- 
ated GW (current epoch) peak frequency is determined 

by cosmology parameters, /max = H*a-ol a * = 1-6 x 
10- 5 Hz (g»/100) 1/6 (T»/100GeV) and is independent of 
turbulence parameters. On the other hand, /max = 
7/max/M is shifted to lower frequency compared to the 
unmagnetized case. From Eq. ([9]), the amplitude of 
MHD-turbulcnce-generated GWs at the peak is a factor 
~ a^ 8 7~ 3 / 4 Af 3 / 4 larger than that in the unmagnetized 
case. 

When modeling turbulence we used the Biskamp and 
Muller model, [30, HH- If we adopt the helical MHD 
turbulence model of Banerjee and Jedamzik [32j (also 
see Refs. [33|) the GW peak frequency remains the same 
while the amplitude of the signal doubles. 

Figure 1 shows that even for small values of magnetic 
helicity the main contribution to the GW energy density 
is from the second, inverse-cascade stage. The GWs will 
be strongly polarized since magnetic helicity is maximal 

. LISA should be able 
2Clj ] . Unlike the unmag- 



11 



at the end of the first stage 
to detect such GW polarization 
netized case due to the second (inverse-cascade) stage 
contribution the GW amplitude is large enough at 10~ 4 
Hz to be detectable by LISA. If the EWPT occurs at 
higher energies (T* > 100 Gev) the peak is shifted to 
higher frequency, closer to LISA sensitivity peak, which 
leads to a stronger signal. Our formalism is applicable 
for GW production at an earlier QCD phase transition, 
assuming the presence of colored magnetic fields 3J| , or 



4 



for any other phase transitions 0]; the peak frequency 
will be shifted according to the changes in T+ and g*. 
The GW signal estimated here exceeds that from bubble 
collisions [i, T^, 3 or from hydrodynamical, unmagne- 
tized turbulence [91 [Tol Il2| . Of course, this strong signal 
assumes initial non-zero (although small) magnetic helic- 
ity so detection of polarized GWs by LISA will indicate 
jarity violation during the EWPT as proposed in Refs. 
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